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INTRODUCTION 

This is the final report, on Grant NGR-34-002-038. It pre- 
sents a review of the activities, of the research carried out 
under the grant, and of the results. It also contains a report 

i 

on the results of research carried out under a no-cost. extension 
during the period July 1, 1968, to March 31, 1969, when the grant 
was terminated. The original grant period was July 1, 1966, to 
June 30, 1968. During this period, the results of the research 
were presented by progress reports. 

One of the objectives of the activities under the grant was 
to assist the National Aeronautics and Space Administration in 
its operations of applying optical methods for data and informa- 
tion handling. This assistance was given by carrying out litera- 
ture searches, by the initiation of short courses on optical data 
processing available to NASA personnel, by carrying out specific 
studies, by research related to pertinent optical processing me- 
thods, such as correlation processing, and by developing new 
methods of analysis and measurement of related optical phenomena. 
The cooperation with the Space Data Control Branch, Code 521 of 
GSFC, and the many stimulating discussions with Mr. Arnold 
Schulman are highly appreciated. 

The research was carried out by the principal investigator 
assisted by graduate and undergraduate students. The names of 
the students working on specific projects were listed in the 
progress reports. One of the students, Dr. J. P. Moffatt, com- 
pleted his dissertation entitled "Effects of Geometric Deviations 
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and Nonlinearities in Coherent Optical Data Processors" in May 
1968 on research under the grant. 
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RESEARCH REVIEW 

The research started with specific problems related to 
the detection of signals in noise by optical correlation pro- 
cessing. Various types of errors which limit the efficiency of 
the optical approach in comparison with electronic methods were 
investigated. The effects of the DC-stop, of finite dimensions 
of the optical signal and filter distributions, and of the non- 
linearities were studied. Attempts were then made to find me- 
thods for the elimination or reduction of these effects. A 
typical example is a function linearizer which is a device to 
be inserted between an electronic signal source and an optical 
signal recorder and which introduces additional nonlinearities 
such that they compensate the original ones. A prototype model 
of such a device was fabricated and operated satisfactorily for 
compensating the nonlinear effects of cathode ray tubes at opti- 
cal signal recording. On the side of this initial research a 
bibliography of books and articles dealing with optical data pro- 
cessing was prepared and published. 

Correlation and spectrum analysis of holograms is another 
subject matter which attracted the interest of students as a 
dissertation topic. Such correlation methods are important since 
they may increase the capabilities of holographic computer mem- 
ories, may be used for three-dimensional-character and body 
recognition, for rendevous-guidance of space vehicles, and for 
the determination of the orientation of various types of bodies. 

A number of elementary steps of the investigation was completed 
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during the grant period. The results were presented in three 
progress reports. 

The study of errors of optical data processors raised the 
question of how to define numerically the quality of the great 
variety of such devices and mfethods. Study of the literature 
and discussions led to adopting a noise- to-signal ratio as a 
characteristic quantity. By considering the deviation of the 
actual optical signal from the hypothetical signal obtained 
under ideal conditions as noise and relating it. to the undis- 
turbed signal intensity, such a quantity, which is customarily 
used in electronic data processing, can also be formulated for 
optical signals and their processing. 

By the use of the noise-to-signal ratio, investigations were 
then carried out to optimize dimensions and operational conditions 
of elements of optical data processors. Tolerances of the dimen- 
sions and geometries of such elements for keeping the quality of 
data processing above specified levels were formulated and found. 
The results of these studies form a part of Mr. Moffatt's dis- 
sertation. 

Considerable efforts were made to find a new method for the 
description of nonlinear effects, since previously used methods 
are inaccurate and difficult to apply as shown in a comparitive 
study. Use of Chebyshev polynomials led to a satisfactory solu- 
tion. The new method is particularly well-suited for optical data 
processing and holography since the signals are mostly sinusoidal 
in these operations or can be represented by such a signal form. 
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An advantageous feature of representing nonlinear char-cteristics 
by Chebyshev series is the fact that expreiments can -be designed 
to check the theoretical results, to find experimentally the 
characteristics of various devices, and to describe them with 
satisfactory accuracy. 

Based on these theories, computer programs were developed 
for the evaluation of the characteristics of photographic plates 
and film, such as spectroscopic plates 649-F, and diagrams made 
which show the effects of nonlinearities of signals recorded on 
such plates. Another program gives the noise-to-signal ratio 
as a function of the bias illumination and the recorded signal 
amplitude. The resulting noise-to-signal ratios presented in 
diagram form indicate the optimum bias levels for obtaining mini- 
mum distortions. The theoretical and computed results were 
checked by experiments and were found in satisfactory agreement. 
The results were presented at meetings of the Optical Society of 
America and found considerable interest indicated by the large 
number of requests for reprints. 

During the no-cost extension, research carried out at a 
reduced level of efforts dealt with a study of nonlinear effects 
on holography. It was found that nonlinearities result in 
secondary or "ghost" images. The theory and an experimental 
verification of these effects is presented in this report. 
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LIST OF INVESTIGATED RESEARCH TOPICS 

1. Preparation of "Bibliography on Optical Information 
and Data Processing" (November 1966) 

2. Effect of the finite width of a DC-stop filter 

3. Nonlinear effects in correlation processing 

4. Methods for the compensation of nonlinear effects in 
optical recording and readout 

5. Effects of the finite width of optical signals and 
filters 

6. Correlation and spectrum analysis of holograms 

7. Effects of the shape of carrier surfaces on optical 
data processing 

8. Unified theory of diffraction 

9. Measurement of the permeability of magnetic fluids 

10. New methods in considering nonlinear effects in optical 
data processing 

11. Representation of optical errors by noise-to-signal 
ratios 

Nonlinear effects in holography 


12 . 
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SEMINARS, CONFERENCES,. AND PUBLICATIONS 

1. Attendance of a short course entitled "Introduction to 
Optical Data Processing," University of Michigan., Ann 
Arbor, Michigan, July 25 to August 5, 1966. 

2. Participation at "Symposium on Modern Optics," Brooklyn 

Inst, of Techn. , New York, March 22-24, 1967. Presented 
paper: "Schlieren Prism." (Co-author, A. T. Shankle) . 

3. Preparation and presentation of a short course on "Mathe- 
matical Fundamentals of Optical Wave Propagation" at 
NASA-GSFC, Greenbelt, Maryland, June 26 to June 30, 1967. 

4. Preparation and presentation of a short course with compu- 
tation sessions on Optical Wave Propagation" NASA-GSFC, 
Greenbelt, Maryland, September 28 to October 4, 1967, 

5. Spring Meeting, Optical Society of America, Washington, 

D. C., March 13-16, 1968. Presented paper: "Nonlinear 

Optical Recording in Data Processing and Holography." 

(J. P. Moffatt and F. J. Tischer) "Unified Theoretical 
Formulation of Electromagnetic Diffraction." (F. J. Tischer) 

6. Publication entitled, "Schlieren Prism for Light Beams" 

A. T. Shankle and F. J. Tischer in "Modern Optics" 
Polytechnic Press, Vol. XVII, 257-264. 

7. Annual Meeting, Optical Soc. of Am. , Pittsburgh, Pa., 

October 9-11, 1968. Presented paper: "Measurement of 

Nonlinear Effects at Optical Recording." 


(F. J. Tischer) 
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8. Spring Meeting# Optical Soc., of Am., San Diego# Calif., 
March 11-14# 1969. Presented paper: "Analysis of Non- 

linear Effects in Holography by Chebyshev Polynomials." 
(P. J. Tischer) 
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LIST OF REPORTS 


The Effect of Finite Width of the DC-Stop in Optical Correlators 
(Preliminary results, 11 pages) August 1966) 

Bibliography on Optical Information and Data Processing 

(Listing of books and articles, 73 pages) November 1966 
The Effect of Non-Linearities - on Optical Correlation Processing 
(Preliminary results, 17 pages) November 1966 
Progress Report (July 1 to December 31, 1966) January 1967 

Subtitles: "The Effects of Finite Spatial Distributions 

and Time Response in Optical Recording." 
"Correlation Analysis of Holograms - General 
Relationships" 

Progress Report (January 1 to April 30, 1967) May 1967 
Subtitles: "Nonlinear Function Modifier" 

"Effects of Shape of Signal Carrier Sur- 
faces in Optical Data Processing" 

"Correlation Processing of Holograms - 
Auto-correlation Processing" 


Progress Report, July 1967 

(Summaries, work in progress, and plans) 

Progress Report (July 1 to September 15, 1967) September 1967 
Subtitles: "Consideration of Nonlinear Effects by 

Chebyshev Polynomials" 

"Expression of the Effects of the Non- 
linearities by a Signal-to-Noise Ratio" 
"Nonlinear Effects in Optical Correlation - 


Continuation" 
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Progress Report (September 15 to December 15, 1967) 30 pages, 
December 15, 1967 

Subtitles: "A Unified Theoretical Formulation of 

Electromagnetic Diffraction" 

"Nonlinear Effects in Optical Processes - 
Continuation" 

"Measurement of the Permeability of a Ferro- 
Magnetic Fluid" 

"Correlation Processing of Holograms - 
Comparison with Spectrum Analysis" 

Progress Report (January 1 to March 31, 1968) March ,1968 
Subtitle: "Final Report on Nonlinear Effects in 

Optical Data Processing" 

Final Report (Period ending June 30, 1968, 161 pages) June 1968 
Subtitle: "Effects of Geometric Deviations and Non- 

linearities in Coherent Optical Data Processors" 
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ANALYSIS OF "GHOST" IMAGES IN HOLOGRAPHY 
BY CHEBYSHEV POLYNOMIALS 

Abstract 

Nonlinear effects in holography resulting in "ghost" 
images at the reconstruction are studied by the use of Cheby- 
shev polynomials. Equations are derived for the transmittance 
of off-axis holograms and for the reconstructed images for the 
two-dimensional case. Using Chebyshev series, the equations 
become relatively simple in the case of line and point objects 
since the various terms of the series represent corresponding 
pairs of the images, of the linear images and of the ghost 
images caused by nonlinearities i Light powers transmitted by 
the various beams can be computed and noise-to-signal ratios 
defined as figures of merit for the holographic processes and 
for the used recording materials. Optimum light levels of the 
reference and signal beams can be found. Experiments are des- 
cribed which confirm the validity of the concepts and of the 
computed relationships. 
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Introduction 

Quality and efficiency of operations in optical data and 
image processing are often degraded by nonlinear effects. Such 
effects are primarily encountered in recording and read-out, but 
also in other operation^ in optical systems. Holography which 
is an optical recording process is a typical example. It is 
strongly affected by nonlinearities. The consequences of the 
nonlinearities are ghost images which are observed at the recon- 
struction. 

The nonlinear relationship in photographic recording can 
be observed at inspecting the characteristic curve for photo- 
graphic material which relates the local amplitude transmittance 
of a, photographic plate or film as the output quantity to the 
exposure as the input quantity. This diagram is more useful from 
a viewpoint of reconstructing holographic images, than the cus- 
tomarily used H & D curve for photographic processing which 
interrelates density and logarithmic exposure. If the former 
relationship were represented by a straight line, higher order 
derivatives being zero, no nonlinear effects would occur and the 
reconstructed image would be identical to the original object, 
disregarding other undesired side effects. This, however, is 
not the case and, depending on the point of operation on the 
characteristic T & E (transmittance versus exposure) curve, a 
varying amount of nonlinear distortions is present which result 
in undesired secondary or "ghost" images at the reconstruction. 
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Nonlinear effects in holography have been investigated by 
1 ” 6 

several authors , however, the derived equations are of con- 
siderable complexity since they usually include multiple corre- 
. 5 

lation integrals - also for simple objects. It then becomes 
difficult to derive a clear understanding of the equations and 
to give a physical interpretation of the relationships between 
the various quantities. It is the purpose of this paper to out- 
line a new method of representation and evaluation of nonlinear 
effects in holography. The method leads to comparatively simple 
relationships particularly at the consideration of rather basic 
object configurations such as lines and points. The resulting 
relationships can be readily interpreted and instructively show 
the origin and optical background of the nonlinear holographic 
phenomena. They also permit finding optimum operational condi- 
tions for keeping these effects small. 

The new method is based on describing nonlinear relation- 
ships between input and output quantities of optical processes 
by Chebyshev polynomials. The figures which describe the 
form of the relationship are the coefficients of a Chebyshev 
series. This approach replaces the usually used methods of re- 
presenting the characteristic relationship by Taylor series, 
Fourier relationships, and gamma-coefficient models. The use 
of Chebyshev polynomials in holography seems to be particu- 
larly advantageous in the case of simple objects. The ghost 
images resulting from nonlinearities are then simple described 
by one Chebyshev coefficient each. In contrast to that, in- 
finite series of coefficients and integrals are required for 
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the description of each image if one of the other known methods 
is used. The described procedure also offers the possibility 
of defining the characteristics of ghost images by individual 
noise-to-signal power ratios and the over-all quality of the 
imaging process by a total noise-to-signal ratio. These quan- 
tities can be computed from the T & E curve of the holographic 
plate o Experiments based on these considerations have confirmed 
the basic concepts and accurate numerical studies of nonlinear 
effects are possible. 
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Method of Analysis 

It is known that an arbitrary continuous relationship be- 
tween two quantities can be described for a finite interval by 
Chebyshev polynomials. A variable u^ is assumed which, between 
its values -1 and +1, represents the finite interval of the input 
quantity. The output quantity u 2 can then be expressed by 

oo 

u 2 “ f(u X> ■ ! b o + L b n W ' (1) 


where the functions T n (u-^) are the Chebyshev polynomials which 
are orthogonal and where the coefficients are given by 


+1 


b - £ /■ 
n 1T _i 


du. 


f(u l> V u l> 


1 - u. 


( 2 ) 


The Chebyshev series has interesting features if the input quan- 
tity varies sinusoidally according to u^ = cos wx. Since the 
terms of the Chebyshev series can be written alternatively as 


one 


T n ( u i) = cos (n cos 1 u x ) 
finds 


(3) 


T n (Uf) = cos nwx , (4) 

and the output quantity becomes as function of the variable x 

u» (x) = rr b + J b cos nwx . (5) 

/ Z o . **, n 

n=l 

The output function is hence a Fourier series with the co- 
efficients equal to those of the Chebyshev series. The first- 
order term with the coefficient b^ has the same form as the 
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input function* The higher-order terms are harmonics caused by 
the nonlinearity of the function f (u^) . 

It is noteworthy that these conditions of a sinusoidally 
varying input quantity are exactly those, encountered in holo- 
graphic processes where interference pattern with such a distri- 
bution of the light intensity are recorded on. the hologram. The 
description by Chebyshev polynomials seems thus ideally suited 
for holographic processing. The input quantity u^ represents 
the varying part of the exposure E which varies sinusoidally 
around a bias value E q as a function of the position on the plate 
described by x for a one-dimensional diffraction pattern. The 
output quantity u 2 stands for the deviation of the transmittance 
from its value for u^ = 0. 

Characteristics of Holographic Plates 
The characteristics of a photographic plate or film for 
holography can be best described by the T & E curves. Examples 
are shown in Fig. 1 for spectroscopic plates 649-F of Kodak, 
which are derived from D and logE curves published by the com- 

7 

pany. The curves can be evaluated by a computer program which 
yields the Chebyshev coefficients according to Eq. (2) . The 
coefficients are computed in terms of two quantities. One is 
the bias exposure E^ which is proportional to the light intensity 
1^ of the reference beam (E^ = I^tjt = exposure time) in the ab- 
sence of a signal beam. The other quantity is the amplitude of 
the sinusoidal variation of the exposure E & due to interference 

with the signal beam. It should be noted that E is smaller 

a 
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than to satisfy the condition that E(x) remains positive. 
The distribution of the exposure is then for a one-dimensional 
interference pattern described by 


E (x) = E, + E cos w x (6) 

D a. X 

The output quantity is the transmittance x(x) = ^(x) which is 
given by the right-hand side of Eq. (5) . The coefficients b R 
are now functions of the bias exposure E^ and the exposure ampli- 
tude E , 
a ' 


b = F (E, ,E ) 
n n b' a 


(7) 


Figure 2 shows as an example a diagram for the first order term 
b 1 which also may be considered as the product of the input am- 
plitude and the average slope of the characteristic curve. In 
the absence of nonlinearities, this term indicates the transfer 
characteristics of the holographic plate. Diagrams similar to 
those shown in Fig. 2 are obtained for the higher-order coeffi- 
cients. If a hologram is illuminated, the higher-order terms 
are related to ghost images where each term of the Chebyshev 
series describes a pair of these, one real and one virtual image. 

Diagrams for the Chebyshev coefficients as shown in Fig. 2 
indicate thus minimum and maximum values of these coefficients 
as a function of the bias exposure and the amplitude of the sinu- 
soidal variation and give an indication of the strength of the 
holographic ghost images. Such diagrams may hence be used for 
the numerical description of the properties of plates and films 
for making holograms. The diagrams can also be used for finding 
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optimum light-intensity levels of the reference beam when dis- 
tortions by ghost images are a minimum. 

An alternative method of describing the quality of holo- 
graphic processes and of plates and films for holography can 
be derived from the concept to consider ghost images as optical 
noise in analogy to similar terms in electronic communications. 
Based on this concept, each ghost image can be described by a 
noise-to-signal power ratio (NSR) which represents the ratio of 
the light power transmitted to form this image divided by the 
power to form the image of the original object. In the case of 
holograms of basic objects such as planes and lines this ratio 
(NSR) simply becomes the ratio of the respective higher-order 
Chebyshev coefficient to the first-order coefficient squared. 

The noise-to-signal ratio of the nth-order ghost is hence 

NSR n = (t^/b ^ 2 

A more general figure is found by adding the relative light 
powers transmitted to form all higher-order ghost images, where 

CO 

NSR total = ( l b n 2 >/ b l * 
n=z 

It seems that this concept facilitates useful comparisons and 
numerical descriptions of processes and materials for holography. 

Theory of Nonlinear Holography 
Based on the preceding equations, relationships can be de- 
rived for the theoretical description of holographic recording. 
The conditions and the geometry which are the basis of the 
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derivation are shown in Fig. 3. A transmitting- type object is 
illuminated by a coherent, quasimonochromatic, and collimated 
light beam. The output beam represents the signal beam which, 
in the plane of the holographic plate, has an amplitude distri- 
bution 

u (x,y) = a(x,y)e j<!,(x ' Y) . (8) 

The reference beam has a distribution at the same location given 
by exp(jw x) , assuming that its transverse amplitude distribu- 

Xi 

tion is constant and normalized to unity. The total exposure 
at the plate is then proportional to the total absolute value 
squared 

2 

E(x,y) = |expjk x x + a (x,y) exp [ j <j> (x,y ) ] | (9) 

Evaluation yields 

E (x,y ) = 1 + a 2 (x,y) + 2a(x,y) cos [k x - <j>(x,y)] / (10) 

V ✓ > ^ 1 \ — ^ 

E b (x,y) E a (x,y) if>(x,y) 

and in simplified writing 

E (x,y) = E b (x,y) + E a (x,y) cos i|^(x,y) . (11) 

The functions E b and E & vary slowly as a function of x and y, 
while cos i p varies fast. The separation of E(x,y) into parts 
varying slowly and fast with position contains an element of 
arbitrariness but it is justified in the case of off-axis holo- 
graphy. In this case the conditions are similar to those in 
electronic modulation theory where similar representations are 

g 

used for modulated- carrier signals. 
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With reference to Eq. (5) , the transmittance resulting from 
the photographic process becomes a linear superposition of terms 

00 

x(x,y) = T Q (x,y) + l x n (x,y) , (12) 

n=l 

where the various terms are given by 

T 0 (-x,y) = x oo [E^ (x, y) ] + j b Q [E b (x,y) , E a (x,y)] , 

T n (x,y) = b n [E b x # y) ,E a (x,y) ] cos mp(x,y) . (13) 

It is noted that the Chebyshev coefficients b n are functions of 
E, and E which in turn are functions of x and y (slowly varying) . 
They represent the envelopes of the rapidly oscillating functions 
cos n^(x,y) and can be obtained from diagrams such as shown in 
Fig. (2) or from the computer program. Reconstruction of the 
first-order term T-^X/y) yields the image of the original object/ 
the reproducable real image and the observable virtual image. 

The higher-order terms of t cause pairs of ghost, images with beams 
linearly superimposed. Since the coefficients b n are functions 
of x and y, they may be represented by Fourier series of the form 

b n (x,y) = l l c cos(ow x x + poo y) . (14) 

o p v y 

It is observed, taking into consideration Eq. (14) , that 
the resulting equations become rather complex in the case of more 
complicated objects, leading to expressions of infinite series of 
multiple integrals. 
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Hologram of a Line Object 

The equations presented in the preceding section become 
more useful if they are applied to simple objects such as a 
line object. The problem becomes then one-dimensional and the 
relationships can be readily interpreted and verified by ex- 
periments. The signal beam originating at the line source 
consists of cylindrical waves which give in the plane of the 
holographic plate a wave amplitude 


W 


1 



exp(jk o 



x 


(15) 


where A^ is proportional to the amplitude of the line source and 
k Q is the free- space wave number or plane-wave propagation con- 
stant. Disregarding the small amplitude variation due the term 
in the denominator , the equation can be written as 


W, = A. exp 


tnw/ 


1 + 


(x/D 1 ) 


- 1 ]} 


(16) 


where the redefined A 1 represents now the amplitude in the sub - 
zenith point of the hologram plane where the origin of the coor- 
dinate system is located as indicated in Fig. 4a. The wave 
amplitude of the reference beam along the plate is 


W 2 = A 2 exp ( jk Q x sinS) , (17) 

under the assumption that the beam hits the plate under an angle 
3 against the plate normal and that A 2 is the amplitude at the 
origin. 

Computation of the intensity of the interference pattern 

* * 

according to I(x) = (W^ + W 2 ) (VI^ + W 2 ) yields 
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I(x) = 2A^ 2 [1 + a 2 + a cos iji(x) ] , (18) 

where 

a = ^ 2^1 ' ~ 2^2 / 

i|)(x) * i k Q _(x 2 /D^) - k x x ? k x = k Q sin 6 

It is noted that Eq. (18) basically is the description of a 
one-dimensional interference pattern with a spatial wave num- 
ber or phase constant in the x direction given by di|>(x)/dx 
which is 

* 

k (x) - difi(x)/dx = k_x/D. " k sin8 . 

X OxO 

Application of Eq. (11) gives the distribution of the exposure 
along the holographic plate 

E (x) = E^ + E(a) cos iJj(x) • (19) 

The bias level which defines the point of operation along the 

2 2 

T & E curve is E^ = A^ t(l + a ) and is constant. The ampli- 
tude of the sinusoidal variation of the exposure E a = A 1 t a 
is constant also. The phase term is given by Eq. (18) 

ip (x) = i k x 2 /D. - k x sin8 • (20) 

2 . O 1 O 

After exposing the plate according to Eq. (20) , and pro- 
cessing it, the distribution of the resulting amplitude trans- 
mittance is given by 


It should be noted that the phase constants in the various 
directions of the local coordinate system are related by 
. 2 , . 2 , . 2 .2 
k x + k y + k z = k o • 
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t(x) = T + l X (x) , 

n=l 


where 


T _ = 


T oo ( V ' + 2 b o (E b' E a> 


( 21 ) 


T 1 = b 1 (E b/ E a ) cos Tjj(x) ; T n = b n (E b ,E a ) cos nij;(x) . (22) 

* 

The Chebyshev coefficients b n are now constants and are inde- 
pendent of the position on the plate. The first order term yields 
the images of the line object if the plate is illuminated by co- 
herent light. The higher-order terms T n (n>l) cause ghost images, 
each term of the Chebyshev series one pair of them. It is ob- 
served that simple relationships are obtained for the line ob- 
ject since each term of the Chebyshev series represents one 
pair of ghost images, one real and one virtual image. This is 
in contrast to other methods of representing nonlinear effects 
such as Taylor series where infinite series are obtained for the 
description of each ghost image. 


General Image Reconstruction 

In considering reconstruction of holograms containing non- 

linearities, Fresnel imaging will be discussed only since 

** 

Frauenhofer imaging seems to be of minor importance. If the 
amplitude distribution of light leaving the hologram is u^(x^,y^). 


5? 

This is an approximation due to disregard of the term in the 
denominator of Eq. (15) . 

* * 

Assuming a hologram size of 2.5 cm squared and a X/4 criterion, 
the distance^ between the image plane and the hologram for Frauen- 
hofer imaging is 4KM (2.5 miles). 
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the distribution in the image plane at a distance D 2 is obtained 

9 

by the modified Fresnel-Kirchhof f diffraction formula in cus- 
tomary notation 


u 2 ( x 2 'Y; 




Ap 


U 1 * x l ,y l* G ^ X 2 


- x l /Y' 


- y 1 )dx 1 dy ] 


(23) 


„ . 4tt q 

K = - j -r— COS 0 , 


G = 


1 exp(-jk o r) 


; r =“^ d 2 2 + (x 2 - x x ) 2 + (y 2 - y-^ 2 , 


4 tt r 

where Q is the angle under which the illuminating light hits the 
hologram and G is the Green's function for free space. 

The light amplitude in the output plane of the hologram is 
proportional to the transmittance and given by 


u i< x i' y i> = A 0 T(x i' y i> = A o > <24) 

n=u 

with t representing the transmittance distribution of the holo- 
gram and where A q is the illuminating light amplitude. The 
various terms of the transmittance are obtained from Eq. (13) , 
where the nth term has the form 

. jni|> (x, ,y 9 ) -jnKx, ,yj 

T n (x l' y l^ = 7 b n^ x l' y l^ + e J (25) 

The two terms between the brackets describe phase variations of 

the hologram which under illumination cause the pair of real and 

virtual images of the nth order. Different terms may cause 

linearly superimposed light beams and images if they hit the 

imaging plane at the same place. The zero-order term represents 

an image of the slowly varying amplitude-envelope distribution of 

the hologram. 
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u 20 {x 2' y 2 } “ K o b o (x l' y l> * < 26 > 

The first-order term describes the original object, 

f f 

u 21 (x 2 ,y 2 ) = K 1 b 1 (x 1 ,y 1 )exp[ji{;(x 1 ,y 1 ) ]G(x 2 -x 1 ,y 2 -y 1 )dx 1 dy 1 

J •> Ap 

(27) 

and the ghost images are given by the higher-order terms 

f 

u 2n (x 2' y 2 )K n b n ^ x l' y l^ exp [ jn^ (x 1 ,y 1 ) ] G (x 2 -x 1 ,y 2 -y 1 ) dx 1 dy 1 

J An 


Reconstruction of the Line Image 
The general relationships of the preceding section are 
applied next to the case of a line object discussed before. 
Under these conditions, the problem becomes one-dimensional and 
the equations can be readily evaluated. Assuming the obliquity 
factor to be unity and disregarding the variations of the dis- 
tance r in the denominator of the Green's function (r - D 2 ) , 

Eq 0 (23) becomes „ 


u 2 (x 2 ) = -j 


-J k o D 2 S -2D 


k x 0 k_ 9 

2D 9 •* 2D 9 ' X 1 zx l x 2 ' 

u, (x, )e dx, . 

Ap X X X 


The amplitude distribution according to Eq. (24) is now given by 


u l (x l> = A o I T n (x l> ' 
n=u 


T n (x l> * b n (E b' E a> te 


jnif) (x.) -jn^U-,) 

+ e ] 


/ 


(31) 
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and 


iKx^) = j k o x l 2//D l - k o x sin ^ * (32) 

For the line object, the Chebyshev coefficients are con- 
stants which depend on the bias exposure E^ and on the amplitude 
of the cylindrical-wave signal beam E as indicated in Eq. (22)* 

a 

The phase term according to Eq. (32) is typical for a, one-dimen- 
sional interference pattern with a linearly increasing phase 
constant (spatial frequency) 


k = = k x./D.-k sing 

x dx o V 1 o 

Considering only the second phase term between the brackets of 
Eq. (31) which yields the real image, one obtains for the; nth 
term of the reconstructed image 


u 2n (x 2> 


Kb 
n n 


k 

o 

2D, 


2 X 1 
(x n -2 x 1 x 2 ) “j k 0 n 2 D" jk Q nx^sing 

e ^ e dx 


Ap 


1 * 
(33) 


Evaluation of the integral for the real image of the original 
object (n = 1) gives under the condition = D 2 = D 


U 21 ^ x 2^ 


= K x b f e 

'Ap 


j (k Q sing- 


x 2 )x l 


dx. 


(34) 


The condition D 2 = implies that the distance between image 
plane and hologram at reconstruction equals that between line 
object and holographic plate when the hologram was made. Further 
evaluation yields 



u 2 i( x 2 ) ~ K<5 (k Q sing - 


(35) 
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The result is a 6-function at the position 

X 2 = D sin 3 (36) 

which represents the image of the original line source. 

Evaluation of Eq. (33) for higher-order values of n shows 
that the ghost images are represented by 6 -functions also, which 
are obtained for positions given by 

D 2 = D-^/n (37) 

and at a transverse distance from the axis 

x 2 = n sin3 = D sin3 . (38) 

The results indicate that the ghost images of the line source 
occur at equal transverse positions x 2 but at different longi- 
tudinal positions Each Chebyshev coefficient and term of 

the series for t describes one of them. 

EXPERIMENTS 

A series of experiments has been carried out to verify the 
derived relationships and to aid at their interpretation. Holo- 
grams were matte of a line object and a point object since the 
described method of analysis is well suited for the study of 
simple objects. The experimental results verify the derived 
equations. 

Experimental Setup 

A major component of the experimental setup, which is illus- 
trated in Fig. 4a, is a HeNe laser with a beam-expanding tele- 
scope which gives a collimated filtered beam. The beam is 
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divided by a beam splitter into two parts of which one passes 
through a cylindrical lens which focuses the beam along a line. 
The beam behind the focal region represents the signal beam 
which has characteristics similar to those of a beam caused by 
a line object. The second part of the beam travels straight 
through the beam splitter and passes through an attenuator. It 
is subsequently reflected by a mirror, and becomes the reference 
beam. Both superimposed beams form an interference pattern which 
is photographed on a Kodak 649-F plate. The photograph repre- 
sents then an off-axis hologram of a, line source. The distance 
between the focus of the cylindrical lens and the plate is 
and the angle under which the reference beam hits the plate is 3. 

Reconstruction is indicated in Fig. 4b. The hologram 
which has a size of about 3/4" by 3/4" is illuminated by the ex- 
panded laser beam. A fraction of the beam passes straight 
through the hologram and forms the zero-order image on the 
photographic plate or film. A second part of the beam is de- 
flected by the holographic interference pattern. of f-axis and 
forms the real image on the plate. A third fraction is de- 
flected into a direction on the other side of the axis and forms 
a beam which gives the impression to originate at the original 
object and thus carries the virtual image. Additional parts of 
the beam are deflected by the higher-order distributions, and 
form the ghost images. Movable masks with rectangular cutouts 
are used to make several holograms and several reconstructed 
images on one plate. 
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Experimental Procedure 

As a first step, a series of holograms was made of a line 
object on 649-F plates at a development time of five minutes 
by the use of the described setup*, The light intensities of 
the reference and signal beams were carefully adjusted to cor- 
respond to a series of operational points along the T&E curve 
shown in Fig. L Six holograms were made on one plate. The 
intensity level of the reference beam was adjusted by inserting 
attenuators of varying attenuation such that six average levels 
of transmittance were obtained ranging from about .8 to .1 
transmittance. The intensity of the signal beam was adjusted 
to correspond to amplitude level of A - 100 (amplitude of sinu- 
soidal variation of exposure) and kept constant when the six 
holograms were made. A. similar set of holograms on one plate 
was made after substitution of a spherical lens for the cylin- 
drical one. The focal region of the lens then represents a point 
source. 

The reconstruction was made by illuminating the holograms 
by an expanded collimated laser beam and by photographing the 
real images on polaroid film. Three pictures were ually taken 
on one frame bymoving an opaque screen with a rectangular cut- 
out up and down in front of the polaroid film holder. The dis- 
tance between hologram and film holder was adjustable to permit 
positioning of the film holder in the planes of the first and 
higher-order images of the line and point object. 
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Results 

The first experiment was made to show first-order and 
higher-order images of a line source. Herewith, one holo- 
gram was used for reconstruction and the real images photo- 
graphed three times in a row with the polaroid film holder 
placed at three distances. The positions coincided with the 
image planes of the original object (line source) and of the 
second-order and third-order ghost images. The three pictures 
are shown in Fig. 5, The upper picture shows a circular illu- 
minated region in the center which represents the zero-order 
circular beam transmitted straight through the circular holo- 
gram. On the right-hand side is the image of the original line 
object. The distance between film and hologram was D 1 which is 
equal to that between object and the holographic plate when the 
hologram was made. The shadow on the left-hand side is caused 
by the beam which, if observed with the eye, would yield the 
virtual image. 

The center picture below shows the second-order ghost 
image of the line object on the right-hand side. The distance 
between film and hologram was reduced to about one-half of 
that when the previous picture was taken. This distance is ob- 
tained by Eq. (37) . It is observed that the lateral distance 
between image and axis is the same as before. This confirms 
the result of the theoretical analysis, Eq. (38). The elipti- 
cal image between the center zero-order image and the second- 
order line image is the result of the first-order beam which, 
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if continuing behind the plane of the film, would focus in the 
image of the original line object photographed in the preceding 
picture above . 

When the lower picture was made, the film holder was moved 
upwards and placed at about D^/3 from the hologram. In this 
plane, the beam resulting from the third-order nonlinearities 
becomes focused and forms a ghost image of the line object. It 
can be noted that the lateral distance from the axis of this 
image is somewhat smaller than in the preceding case. The 
elliptical images at the left of the line image are the result 
of the first and second order beams. It is observed that the 
light intensities of the reconstructed line images considerably 
decrease with increasing order number. The corresponding light 
amplitudes are proportional to the individual Chebyshev coeffi- 
cients for the specific light levels when the hologram was made. 

The purpose of the second experiment was the study of the 
influence of the bias level of the reference beam on the non- 
linearities. The intensities of the second-order ghost images 
were compared at various amplitude levels of the reference 
beam. The amplitude of the signal beam was kept constant as 
indicated in the description of the experimental procedure in 
section B. The hologram was that of a line object. Figure 6 
shows the results. Photographs were made on polaroid film of 
the real images of three holograms by moving the film holder 
up and down. The pictures were taken with the film plane in 
the image plane for the second-order ghost images. The ghost 
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images are focused in this plane and have the form of the origi- 
nal line object. The line images are shown on the right-hand 
side. The beams for the linear real and virtual images of the 
line object cause the elliptical images on both sides of the 
centrally located zero-order image. 

In the upper picture, the light level of the reference beam 
was in the region for minimum noise-to-signal ratio (NSR) causing 
a transmittance of about ,5 and giving minimum nonlinearities. 

The ghost image is rather weak. The lower picture shows the re- 
constructed images when the level was adjusted for a transmittance 
of about .15 giving a maximum NSR and hence large nonlinear effects. 
Correspondingly the intensity of the ghost image is rather strong. 
The picture in the center corresponds to an operational point on 
the T&E curve between the maximum and minimum at a transmittance 
of about .3 (see Fig. 1). At the reconstruction, the light level 
of the first-order beam was kept constant. 

In the third experiment, the cylindrical lens was replaced 
by a spherical lens which formed a point source in the focal 
region. The resulting interference pattern is a hologram of a 
point object. Reconstruction in a slightly expanding beam re- 
sulted in the three photographs shown in Fig. 7. The upper photo- 
graph shows the reconstructed real image in the plane of the 
linear image where the right-hand illuminated point represents 
the image of the original object. The circular illuminated region 
in the center is the image resulting from the zero-order, straight- 
through beam. The two photographs below show the ghost images 
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caused by the second-order and third-order nonlinearities. The 
circular images between the point images on the right-hand side 
and the zero-order images in the center are the result of illumi- 
nation by the lower-order beams. Light shadows left of center 
are the result of beams carrying the virtual image. 

Conclusions 

It was shown that relationships can be derived by the use 
of Chebyshev polynomials which describe the nonlinear effects in 
holographic processes. Equations were found for the transmit- 
tance of thin off-axis holograms of general, two-dimensional ob- 
jects and for the reconstructed images in the Fresnel zone. The 
equations are complicated for general objects, as is the case at 
using other methods of representing nonlinear characteristics, 
such as Taylor series. The equations, however, become manageable 
for fundamental objects, particularly for line objects where the 
problem is one-dimensional. The results can be readily interpreted 
and were confirmed by experiments. 

The equations and corresponding experiments confirm the pre- 
viously reported observations that the nonlinear effects cause 
secondary or "ghost" images. The use of Chebyshev polynomials 
for the description of nonlinear characteristics offers the ad- 
vantage that in the case of fundamental objects, the description 
of the transmittance distribution of the hologram becomes rather 
simple. Each term of the Chebyshev series directly represents a 
pair of ghost images. This is not the case if other methods of 
describing nonlinearities are used. Infinite series are then 
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obtained instead of single terms „ If the hologram is illumi- 
nated by coherent lights the elementary harmonic distributions 
which superimposed represent the transmittance and which are 
described by the various Chebyshev coefficients deflect fractions 
of the beam passing through the hologram,, The fractional beams 
are focused at different positions and form the images . The de- 
scribed experiments follow instructively the derived equations 
as shown by the reproduced photographs of the real image in 
Figs 0 5 through 

Another merit of the new method is the fact that the non- 
linear effects can be treated numerically „ Using Chebyshev 
polynomials , the various coefficients of the series represent 
pairs of ghost images * Hence, the light power carried in the 
beams forming the various images actually can be computed 0 
Ratios for the light powers forming the various ghost images 
to that of the actual linear image can be foundo This permits 
determination of individual noise-to-signal ratios and a cor- 
responding over-all ratio „ The latter represents the ratio of, 
the light power forming all ghost images to that of the linear 
image of the original object „ Hence, one can find figures of 
merit for the quality of holographic reproduction with regard 
to nonlinearities o This in turn opens the way for the numeri- 
cal evaluation and comparison of film, plate, and other recording 
material for holography and of the related processing methods „ 

It also permits determination of optimum operational conditions 
given by the ratio of the level of the signal beam to that of 
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the reference beam. The new method thus seems to have a 
of merits, however, further studies and further use will 
an indication of their importance. 


number 

give 
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Characteristic T&E curve (transmittance-versus-exposure) for 649-F Spectroscopic plates 
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